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In the recent years, numerous studies have been done on Stirling cycle and Stirling engine which have 
been resulted in different output power and engine thermal efficiency analyses. Finite speed thermo¬ 
dynamic analysis is one of the most prominent ways which considers external irreversibilities. In the 
present study, output power and engine thermal efficiency are optimized and total pressure losses are 
minimized using NSGA algorithm and finite speed thermodynamic analysis. The results are successfully 
verified against experimental data. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

One of the greatest issues for the societies in the near future is 
security of sustainable energy [1]. Great amount of enterprise will 
be vital to satisfy future energy requirements with climate-friendly 
energy technologies [2-5], comprises electricity produced from 
solar and wind energy [6-8], bio-fuels [9,10], and carbon dioxide 
(CO 2 ) capture and storage (CCS) [11—14]. In this regard, scrutinized 
planning with respect to the prospect plays an important role to¬ 
ward better future. 

Carbon dioxide emission is one of the restrictions of applying 
fossil fuels in all over the world which motivates expansion of 
higher efficiency and maximal electric energy per 1 kg CO 2 emitted 
technologies [15,16]. Use of Stirling cycles can lead to reduce the 
CO 2 emissions to the atmosphere. 
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Stirling cycle is one of the main primitive standard air cycles for 
heat engines [17,18]. External heat source and high efficiency are 
the Stirling engines advantages. Stirling engines can use solar en¬ 
ergy which is available in one-third of the day, in the result solar/ 
fuel hybrids are recommended. The combustion of the Stirling 
engine is continuous process and can burn all types of fuel with any 
quality [18-20]. The Stirling engine can theoretically be a highly 
efficient engine to convert heat into mechanical work at the Carnot 
efficiency when the ideal regeneration, isothermal compression 
and expansion are considered. The thermal limit for the operation 
of a Stirling engine depends on the working temperatures of the 
heater and cooler sides. In most instances, the engine operates with 
a heater and cooler temperature of 923 and 338 K, respectively [21 ]. 
The engine efficiency varies from about 30 to 40% which results in 
forming a typical temperature range of 923-1073 K, and normal 
operating speed range from 2000 to 4000 rpm [22-27]. The in¬ 
fluence of the regenerator effectiveness, the dead volume on the 
power output and the thermal efficiency are studied by Kong- 
tragool and Wongwises [24]. 

Carlson et al. provided a non-isothermal heat exchanger model 
which eliminates the necessity for infinite heat transfer time 
correspond to slow engine speed that is impractical [28]. 
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Organ studied the effects of various parameters such as diam¬ 
eter, length and materials on regenerator performance, irrevers¬ 
ibilities and temperature gradient in Stirling engine regenerator 
while the regenerator is optimized [29,30]. 

Martaj et al. presented a thermodynamic analysis of a low 
temperature differential Stirling engine at steady state operation, 
and energy, entropy and exergy balances were presented at each 
main element of the engine [31]. 

Formosa and Despesse modeled engine power output and ef¬ 
ficiency due to dead volume by implementing of the isotherm 
model [32]. 

Timoumi et al. developed a numerical model based on lumped 
analysis approach which takes into account losses and is imple¬ 
mented for optimization of GPU-3 Stirling engine [33-35]. 

The effects of heat transfer, regeneration time, and imperfect 
regeneration on the performance of the irreversible Stirling engine 
cycle are investigated by Wu et al. [36]. Li et al. developed a 
mathematical model for the overall thermal efficiency of solar 
powered high temperature differential dish Stirling engine with 
finite heat transfer and the irreversibility of regenerator and opti¬ 
mized the absorber temperature and corresponding thermal effi¬ 
ciency [37]. Tlili studied the effects of regenerating effectiveness 
and heat capacitance rate of external fluids in heat source/sink at 
maximum power and efficiency [38]. The effects of addressed ir¬ 
reversibilities of regeneration and heat transfer of heat/sink sources 
were investigated by Kaushik and Kumar [39,40]. 

This paper concentrates on the evaluation of the losses because 
of the irreversibilities in the Stirling cycle which is an issue of 
substantial interest to those concerned with the performance and 
analysis of heat engines. Our researches and latest studies by 
others [41,42] have scrutinized that irreversibilities in the ther¬ 
modynamic cycle have a specific significance in forecasting the 
performance of Stirling engines. Remarkable attempts have been 
made over the past years to completely recognize the irrevers¬ 
ibilities associated losses [43,44]. Finite Speed of the piston, Fric¬ 
tion and Throttling in the regenerator cause the internal 
irreversibilities which have to be taken into account for validating 
of any design of solar Stirling engines computation. One of the 
earliest efforts for a validation of the solar Stirling engines scheme 
of computation has been done by Costea et al. [45]. More exactly 
validation has been possible to accomplish through a validation of 
the schemes of computation for as many as possible Stirling en¬ 
gines. Probably, it can be done by using the new branch of irre¬ 
versible thermodynamics, so-called “Thermodynamics with Finite 
Speed, and the Direct Method” [46]. 

Based on this, in 2002 it became possible to obtain the most 
powerful validation of a scheme of computation of the perfor¬ 
mances (efficiency and power) for 12 Stirling engines (the most 
performing in the world) working in 16 regimes [47]. Only with 
such a powerful “tool” it is possible to obtain a validation for a 
scheme of computation of solar Stirling engines performances for 5 
solar Stirling engines [48-50]. 

This method has been employed in a number of models for the 
analysis and optimization of Stirling engines including the impact 
of irreversibilities in the cycle, successfully. Flowever, analyses 
that are based on the entropy generation or exergy techniques do 
not relate the irreversibilities to the physical phenomena that 
cause them. The model presented here directly associates the ir¬ 
reversibilities to the operation of the cycle at finite speed. On the 
basis of the entropy generation techniques, for the studying of 
Stirling engine cycle performance, Costea et al. [48,49] have 
included the impacts of heat transfers, imperfect heat regenera¬ 
tion and irreversibilities of the cycle including pressure loss 
associated with the finite speed of the piston and Displacer 
throughout processes as gas passes through the regenerator, 


heater and cooler, and finally mechanical friction due to the mo¬ 
tion of the moving parts. 

For considering all mentioned issues, we have studied three 
objective functions: output power, the thermal efficiency of the 
entire solar Stirling system and pressure losses. In Addition, the 
multi-objective optimization is conducted with eleven decision 
variables including the temperatures of heat source and heat sink 
and their difference with working fluids, rotation speed, mean 
effective pressure and stroke. 

Solution of the multi-objective optimization problems is an 
extremely difficult goal which requires the simultaneous satisfac¬ 
tion of a number of different and even conflicting objectives. 
Evolutionary algorithms (EA) were initially extended and employed 
during the mid-eighties in an attempt to stochastically solve 
problems of this generic class [51 ]. A reasonable solution to a multi¬ 
objective problem is to investigate a set of solutions, each of which 
satisfies the objectives at an acceptable level without being domi¬ 
nated by any other solution [52]. Multi-objective optimization 
problems, in general, show a possibly uncountable set of solutions 
namely as Pareto frontier, whose evaluated vectors represent the 
best possible trade-offs in the objective function space. In this term, 
multi-objective optimization of different thermodynamic and en¬ 
ergy systems has been paid attention by researchers nowadays 
[53-58]. 

2. Stirling system 

As shown in Fig. 1, Stirling cycle consists of four major processes. 
Process 1-2 is an isothermal process, in which the compressing 
working fluid rejects the heat at constant temperature (Tc) to heat 
sink which has a constant temperature (Tl). Then the working fluid 
crosses over the regenerator and is warmed up to Th in an isochoric 
process 2-3. In process 3-4, the working fluid expands at a constant 
temperature, Th, and obtains the heat from the heat source at a 
constant temperature ( T H ). Last process (4-1), is an isochoric cool¬ 
ing process, where the regenerator absorbs heat from the working 
fluid. In an actual cycle it is impractical to have an ideal heat 
transfer in the regenerator in which the entire amount of absorbed 
heat (in the process 4-1) is transferred to the working fluid into the 
isochoric heating process (process 2-3). 

3. Analysis of the Stirling engine cycle with irreversibilities 

The pressure losses presented in this task is as follows [45-49]: 

A Pi = A Pthrot R + A Pf + A Pw ( 1 ) 

where Ap throtR is the pressure drop coming of the internal friction of 
the current which takes place in the regenerator and is insignificant 
in heat exchangers (coolers and heaters) [45-50]. 



Fig. 1 . P—V and T—S diagrams of an isothermal Stirling engine cycle [48]. 
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Ap/is the pressure drop which comes of the mechanical friction 
of the parts of the engine and has been experimentally [45-49] 
obtained: 
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And finally the pressure drop coming of the piston speed 
[45-49] is calculated as follows: 
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A Qr is the heat loss during the two regenerative processes in the 
cycle which can be obtained from Eq. (2) [45-47]: 


A Qr = m g C Vg X(T H - A T H — T L — A T L ) 


( 2 ) 


The heat released between heat source and working fluid (Qh) is 
obtained as follows [45-47]: 


Qj = m (A — + l)(i + 1) _ bAp throtR _ fkp/ \ 

^ ~ g V 4p m 2 p m p m J 

R(T h - AT H )ln X 


( 3 ) 


The net heat released from the heat source (Qh) is obtained as 
follows [45-47]: 


Qh = Qh + A Qr 


( 4 ) 


The pure flux transferred to working fluid is obtained as follows 
[45—47]: 


Qh = (Qh + ^Qr) 
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Substituting Eqs. (2)—(4) into Eq. (5), the pure flux transferred to 
working fluid, we have [46-52]: 
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The output power of the engine is calculated as follows [45]: 
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By placing Eqs. (7a)—(7c) in Eq. (7), we have: 


(7c) 
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(8b) 


where p c is Carnot efficiency and r]\\,irr is the efficiency of the second 
law. In this evaluation, each one of the irreversibilities has an in¬ 
dependent effect which is taken into consideration in the efficiency 
of the second law, in such a way that [45-50]: 


Vll.irr — VlI,irr(X )' VlIJrr(Ap) 


( 9 ) 


The effect of imperfect regenerating is calculated in TjiiMX) 
which the following are [45-49]: 
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where y is the second adjusting coefficient of the method, and Xi, X 2 
are regenerative losses coefficient values corresponding to an 
“optimistic”, respectively “pessimistic” evaluation of them. The 
value of y has been determined from experimental data of 4 solar 
Stirling engines [47]. Settingy = 0.72 provided the best fit between 
analytical results and experimental one, so this value is carried out 
as an adjusting coefficient in the theoretical analysis. In the above 
relation [46-49]: 
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In which the convection coefficient h is computed as [46-49]: 
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On the other hands [47]: 
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The effect of mechanical friction, the speed of the piston and the 
pressure drop in regenerator are calculated in ^ii )ir r(Ap) [45-49]: 
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The above analysis shows that the pressure losses and their 
effect on efficiency and output power of the engine depend on the 
mean piston speed and hence the engine rotation speed. 


4. Multi-objective optimization with evolutionary algorithms 



4.1. General concepts of the multi-objective optimization 


Fig. 2. Schematic of the objectives space. 


Consider a decision-maker who wants to optimize I< objectives 
that objectives are non-commensurable and the decision-maker 
has no clear priority of each objective over the others. Without 
loss of generality, all objectives are of the minimization type. For a 
minimization type objective, it can be renewed to the maximiza¬ 
tion type by multiplying negative one. A minimization multi¬ 
objective decision problem with I< objectives is defined as follows: 

Given a non-dimensional decision variable vector x = {xi,...,x n } in 
the solution space X find a vector x* that minimizes a given set of K 
objective functions f[x) = {/i(x*),...,//<(x*)}- The solution space X is 
generally restricted by a series of constraints, such as gfx*) = bj for 
j = l,...,m, and bounds on the decision variables [52]. Overall, no 
solution vector X exists that minimizes all the K objective functions 
simultaneously. In another word, in various real-life problems, 
objectives at issues conflict each other. So optimizing X with respect 
to a uni-objective often leads to unacceptable results with respect 
to the other objectives. Thus a new concept, well-known as the 
“Pareto optimum solution”, is implemented in multi-objective 
optimization problems. A feasible solution X is called “Pareto 
optimal” if there is no other feasible solution Y that dominates 
solution X. By definition, a feasible solution Y is said to dominate 
another feasible solution X, if and only if, f{Y) < f(X) for i = 1,...,/< 
and/,( Y) < fj(X) for at least one objective function. This means that a 
feasible vector X is called Pareto optimal if there is no other feasible 
solution Y that would reduce some objective function without 
causing a simultaneous increase in at least one other objective 
function. The set of all feasible non-dominated solutions in X is 
referred to as the “Pareto optimal set”, and for a given Pareto 
optimal set, the corresponding objective function values in the 
objective space are called the “Pareto optimal frontier” [52]. 

Notice Fig. 2 for having an effective insight into the concept of 
Pareto optimal. The objective function space is illustrated for an 
optimization problem with two objective functions f\ and f 2 . The 
feasible and infeasible areas are shown as the areas inside and 
outside of the curve respectively. Each point in the feasible area can 
be represented as a solution of the problem. In Fig. 2, the values of 
both functions f\ and f 2 . For point M is lower than the corre¬ 
sponding values of the point j. Thus the point M dominates point J 
or point M is better than a point J. In the same way points L and N 
dominate M. But points like I and K neither dominate M nor are 
dominated by M. Therefore, only those points that are located in the 
left-down parts of M will dominate M (like L and N). Now consider 
point R, there is no point in the feasible area located in the left- 


down part of it. Thus R is not dominated by any other feasible 
point. Therefore R is a Pareto optimal. This is true for points P, A, R , E, 
T, 0 and all of the other points located at the bold curve indicated as 
“The Pareto Optimal Frontier” in Fig. 2. Notice that, in the feasible 
area, the minimum values of /j and f 2 belong to points P and 0 
respectively. Thus P and 0 are the solutions for uni-objective opti¬ 
mization problems which their objective functions are f\ and f 2 , 
respectively. Other points on the Pareto frontier are also optimal, 
but each of them should be chosen as the final solution which is 
required a decision-making process. 

4.2. Optimization via EA 

In the present study, the Pareto frontier is achieved using the 
Genetic Algorithm (GA) known as a tributary of evolutionary al¬ 
gorithm. Genetic Algorithms were evolved by John Holland in the 
1960s for involving the mechanisms of natural adaptation into 
numerical optimization and computer algorithms [54]. They are 
employed as a computer simulation in which a population of ab¬ 
stract representations (called chromosomes or the genotype of the 
genome) of candidate solutions (called individuals, creatures, or 
pheno-types) to an optimization problem develops toward better 
solutions. The progress often begins from a population of randomly 
generated individuals and happens in generations which the fitness 
of every individual in the population is evaluated In each genera¬ 
tion; multiple individuals are stochastically selected from the 
existing population (based on the fitness), and modified (recom¬ 
bined and possibly randomly changed) to create a new population. 
The new population is applied in the next iteration. Generally, the 
algorithm ends when either a satisfactory fitness level has been 
overtaken for the population, or a maximum number of genera¬ 
tions have been made. If the algorithm has finished owing to a 
maximum number of generations, a satisfactory solution may or 
may not be found. In genetic algorithms, a candidate solution to an 
issue is routinely entitled a chromosome, and the evolutionary 
livability of every chromosome is represented by a fitness function 
which is a powerful optimization method for nonlinear problems 
[53,55]. 

Furthermore, multi-objective evolutionary algorithms (MOEAs) 
have been expanded over the recent decades by many tests on 
twisted mathematical problems and on real-world engineering 
subjects, they can remove the obstacles of classical procedures 
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[53,55]. The structure of the MOEA employed in the present 
investigation by the authors is shown in Fig. 3 [54]. The real 
values of decision variables are taken into account instead of their 
binary coded. 


0.05 < D c < 0.14 m 

(24) 

0.02 <D r < 0.06 m 

(25) 


43. Objective functions, decision variables and constraints 


5. Decision-making in the multi-objective optimization 


Three objective functions for this study are the system pressure 
losses (]TAPj), the output power (Power) and the Stirling engine 
thermal efficiency ( 77 ), respectively. 

In this paper eleven decision variables have been considered as 
follows: 

n r : Engine’s rotation speed 
Pm- Mean effective pressure 
s: Stroke 

Th‘. Temperature of heat source 
T l : Temperature of heat sink 

A T l : Temperature difference between heat source and working 
fluid 

A Th'. Temperature difference between heat sink and working 
fluid 

Nr: Number of gauzes of the matrix 
Dc‘. Piston Diameter 
Dr: Regenerator Diameter 
L: Regenerator’s length 

The objective functions with respect to following constraints 


have been solved: 


1200 <n r < 3000 rpm 

(15) 

0.69 < P m < 6.89 MPa 

(16) 

0.06 < s < 0.1 m 

(17) 

800 < T h < 1300 I< 

(18) 

288 < T l < 360 K 

(19) 

64.2 < A T h < 237.6 K 

(20) 

5 < AF l < 25 K 

(21) 

250 < N r < 400 

(22) 

6 x 1CT 3 < L < 73 x 1CT 3 m 

(23) 



Fig. 3. Scheme for the multi-objective evolutionary algorithm used in the present 
study. 


A procedure of decision-making is required in multi-objective 
optimization to choose the final optimal solution from available 
solutions. Great numbers of techniques for decision-making pro¬ 
cess in decision problems exist. The methods can be utilized in 
decision-making to opt a final optimal solution from the Pareto 
frontier which is gained through optimization. Insomuch dimen¬ 
sion of various objectives in a multi-objective optimization issue 
might be various (for example in this work the thermal efficiency 
objective has no dimension while the dimension of the output 
power is in Watt). It is true that, dimensions and scales of objectives 
space have to be unified at first. In turn, objectives vectors have to 
be non-dimensionalized before the decision-making procedure. 
There are some ways of non-dimensionalization utilized in deci¬ 
sion-making including linear, Euclidian, and fuzzy. 

• Linear non-dimensionalization 


Consider the matrix of objectives at various points of the Pareto 
frontier is denoted by Fy where i is index for each point on the 
Pareto frontier and j is the index for each objective in objectives 
space. Therefore a non-dimensionalized objective, Fy, is defined as, 




max (. Fy) 


For maximizing objective 


(26a) 



max (1 /Fy) 


For minimalizing objective 


(26b) 


• Euclidian non-dimensionalization 


In this method, a non-dimensionalized objective, Ffj, is defined as, 
F • 

Fjj = . IJ For minimizing and maximizing objectives 

yjYXL 1 (fy) 2 

(27) 

• Fuzzy non-dimensionalization 


In this method, a non-dimensionalized objective, Ffj, is defined 
as [55,56], 



Fy ~ min (Fy) 
max (Fy) - min (Fy) 


For maximizing objectives 


(28a) 



max (Fy) - Fy 
max (Fy) - min (Fy) 


For minimizing objectives 


(28b) 


In this work, most recognized and common types of decision¬ 
making processes including the fuzzy Bellman-Zadeh, LINMAP 
and TOPSIS methods are implemented in parallel and the final 
optimal solution has been decided based on engineering experi¬ 
ence and criteria among solutions which is suggested by these 
three methods. The Bellman-Zadeh method utilizes the fuzzy non- 
dimensionalization while the other methods (LINMAP and TOPSIS) 
employ Euclidian non-dimensionalization. The following sections 
are presented here in order to describe these decision-making 
algorithms. 
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5 A. Bellman-Zadeh decision-making method 


For maximized objective functions: 


When using the Bellman-Zadeh approach, each F/X) is replaced 
by a fuzzy objective function or a fuzzy set, 

Aj = {X, fi Aj (X)} XeL, j = 1,2, ...k (29) 

where g Aj {X) is a membership function of Aj [59]. 

A final decision is defined by the Bellman and Zadeh model as 
the intersection of all fuzzy criteria and constraints and is repre¬ 
sented by its membership function. A fuzzy solution D with setting 
up the fuzzy sets (Eq. (30)) is turned out as a result of the inter¬ 
section D = nf =1 Aj with a membership function, 


MFiPO — < 


1 if Fj(x) > F™ ax , 

Fj - F™ in 


cmax _ cmin 
i i 


0 


if F™ in < Ff < F™ ax , 
if F/ (x) < F™ in 


Fuzzy constraints are formulated: 


GjWKGP^ + dj, j = 1,2, 


(35b) 


(36) 


MdW = H f^n A (X) = min u A (x), XeL (30) 

Using Eq. (28a), it is possible to obtain the solution proving the 
maximum degree as follows, 


where dj is a subjective parameter that denotes a distance of ad¬ 
missible displacement for the bound Gf 1 ** of the jth constraint. 
Corresponding membership functions are defined in following 
manner [55,56]: 


ma xfi D (X) = max min n A (x) 

XeL J = l, J 

(31) 

X° = arg max min g A . (x) 

XeL j= !,•••,/< J 

(32) 


To obtain Eq. (31), it is necessary to build membership func¬ 
tions gAj(X ); j = 1, reflecting a degree of achieving “own” 
optimal by the corresponding Fj(X), X c L; j = 1 ,...,/<. This 
is satisfied by the implement of the membership functions [60]. 
The membership function of objectives and constraints, linear 
or nonlinear, can be chosen depending on the context of problem. 
One of possible fuzzy convolution schemes is presented 
below [61]. 

Initial approximation for X-vector is chosen. Maximum (mini¬ 
mum) values for each criterion Fj(X) are established via scalar 
maximization (minimization). Results are denoted as “ideal” points 
{X?, j = 1.m}. 

The matrix table {T}, where the diagonal elements are “ideal” 
points, is defined as follows: 



F 2 (x?)... F„(x°)i 
F 2 (x°) ••• fn(x 2 °) 


(33) 


G/(x) > G™ ax 

G™ ax < G/(X) < G™ ax + dj 
G/(x) < Gp ax 

(37) 

A final decision is determined as the intersection of all fuzzy 
criteria and constraints represented by its membership func¬ 
tions. This problem is reduced to the standard nonlinear pro¬ 
gramming problems: to find such values of X and k that 
maximize k subject to 


Fci(X) 


0 

Gj(x) - C™x 

1- 


dj 


1 


if 

if 

if 


A < /ip., i = 1,2,..., n 
^ < FGj» j = 1,2, ...,/c 


(38) 


The solution of the multi-criteria problem discloses the meaning 
of the optimality operator and depends on the decision-makers 
experience and problem understanding. 


5.2. LINMAP decision-making method 


LFi(x°) F 2 (x„)... F„(x°)J 

Maximum and minimum bounds for the criteria are defined: 
pmm = m j n p, ^ I = \ n (34a) 

F™ = max Fj (x ^, i = 1,..., n (34b) 


An ideal point on the Pareto frontier is the point in which 
each objective is optimized regardless to satisfaction of other 
objectives. It is clear that in the multi-objective optimization it 
is impossible to have each objective in its optimal condition 
which can be obtained only in a single-objective optimization. 
Therefore, the ideal point is not located on the Pareto frontier. 
In the LINMAP method, after Euclidian non-dimensionalization 
of all objectives, the distance of each solution on the Pareto 
frontier from the ideal point denoted by d 1+ is determined as 
follow, 


The membership functions are assumed for all fuzzy goals as 
follows. 

For minimized objective functions: 




pldeal 

j 


(39) 




0 


F™ x - Fj 

cmax _ cmin 
i i 


if Fj (x) > F^ ax , 


if F™ 1 < Fj < F^ ax , 


1 


if Fj(x) < F™ n 


where n denotes the number of objective while stands for each 
solution on the Pareto frontier (i = 1,2,...,m). In Eq. (39), Fj deal is the 
ideal value for jth objective obtained in a single-objective optimi¬ 
zation. In LINMAP method, the solution with minimum distance 
from ideal point is selected as a final desired optimal solution, 
hence, i index for a final solution, ifi na i is, 
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ifinai = is min (dj + ); i = 1,2(40) 

5.3. TOPSIS decision-making method 

In this method beside the ideal point a non-ideal point is 
defined. The non-ideal point is the ordinate in objectives space in 
which each objective has its worst value. Therefore, beside the 
solution distance from ideal point, d I+ , the solution distance from 
non-ideal point denoted by d z -_ is utilized as a criterion for selection 
of the final solution. Hence [55,56], 




pnon 

j 



(41) 




Fig. 5. Pareto optimal frontier in objectives’ space (77 - Power). 


In continuing the TOPSIS method a Clj parameter is defined as 
follows, 



dj- 

di+ + d\_ 


(42) 


In TOPSIS method a solution with minimum C/ z is selected as a 
desired final solution, therefore, if ip na i is index for the final selected 
solution, we have, 


x 10 



Fig. 4. (a) Pareto optimal frontier in objectives’ space, (b) Pareto optimal frontier in 
objectives’ space. 


Fig. 6. Pareto optimal frontier in objectives’ space (^Ap,- - 77 ). 
















































320 


M.H. Ahmadi et al. / Renewable Energy 60 (2013) 313—322 


x 10 



Fig. 7. Pareto optimal frontier in objectives’ space (^Ap, - Power). 


ifmai = max (C/j); i = 1,2, ...,m (43) 


6. Results and discussion 

The output power and thermal efficiency of the Stirling engine 
are maximized simultaneously and the pressure losses of the Stir¬ 
ling engine are minimized simultaneously using the multi¬ 
objective optimization method which works based on the NSGAII 
algorithm. 

In this regard, optimization is performed with objective func¬ 
tions that are expressed by Eqs. (1), (7) and (8a) constraints which 
are expressed with Eqs. (15)—(25). 

Design parameters (decision variables) of optimization are the 
engine’s rotation speed, Mean effective pressure, Stroke, Number of 
gauzes of the matrix, Piston Diameter, Regenerator Diameter, 
regenerator length, temperature of the heat source, temperature of 
the heat sink, Temperature difference between the heat source and 
working fluid, Temperature difference between heat sink and 
working fluid. 

In order to have consistency with previous works, specifications 
of the Stirling engine are considered as follows [35]: 

N = 8, v = 3.249 x 1(T 5 (m 2 s” 1 ), y = 1.667, C v = 3115.6 
0 kg” 1 [<-’), Pr = 0.71, m g = 0.001135 (kg), p st = 8030 (kg nr 3 ), 
C Ps = 5193 0 kg- 1 K- 1 ), C R = 502.48 (J kg- 1 K _1 ), A = 1.2, 
b = 6.88 x 10- 5 ,/= 0.556, d = 4 x 10“ 5 (m). 


Table 2 

Error analysis based on the mean absolute percent error (MAPE) method. 


Decision-making 

Bellman—Zadeh 

LINMAP 


TOPSIS 



method 










Objectives 

Power 

VSE 

EA Pi 

P 

VSE 

EA Pf 

Power 

VSE 

EA Pi 

Average error (%) 

2.70 

2.90 

5.30 

2.80 

3.80 

4.60 

3.00 

2.00 

6.30 

Max error (%) 

3.80 

4.60 

10.50 

5.00 

6.60 

10.40 

8.50 

2.60 

10.70 


Pareto optimal frontier for three objective functions, output 
power, thermal efficiency of Stirling engines and system pressure 
losses are represented in Fig. 4a. The chosen points based on de¬ 
cision-making methods are shown, too. The equation of the Pareto 
optimal frontier is obtained and indicated in Fig. 4b, which can be 
applied to gain system pressure losses, as following: 

AP,- = 1.736 x 10 4 + 7795p+4517P +1597r/P+3472p 2 +77.77P 2 
R 2 = 98.45 

(44) 

Pareto optimal frontier for output power and thermal efficiency 
are shown in Fig. 5, which are two objective functions. It is clear 
that output power is between 5300 W and 8300 W in the curve 
fitted diagram: 

rj = -2.23 x 10" 14 P 3 - 1.494 x l(r 10 P 2 - 3.917 x 10“ 6 P 
+ 0.204 

(45) 

Pareto optimal frontier for system pressure losses and efficiency 
are shown in Fig. 6, which are two objective functions. It is obvious 
that thermal efficiency is between 0.03 and 0.18 in the curve fitted 
diagram: 

A Pi = 9294e 1497 ' 1 + 0.006339e 7915 '' (46) 

Pareto optimal frontier for system pressure losses and output 
power are shown in Fig. 7, which are two objective functions. It can 
be seen that output power is between 500 W and 8500 W in the 
curve fitted diagram: 

AP, = 1.991 x 1(T 10 P 3 6 + 1.026 x 10 4 (47) 

The optimal results for decision parameters and objective 
functions using LINMAP, TOPSIS, Bellman-Zadeh decision-making 
methods are summarized in Table 1. 

Maximum and mean error analyses for objective functions are 
categorized in Table 2. The analysis is done for all mentioned 
decision-making methods. 

Sensitivity analysis for volumetric ratio and wire diameter are 
shown in Tables 3a and 3b based on the highest and lowest values. 
Decision variables are represented in Table 3a using TOPSIS 
decision-making method. The results of sensitivity analysis on 
objective functions are indicated in Table 3b using TOPSIS decision¬ 
making method. 


Table 1 

Decision-making of multi-objective optimal solutions. 


Optimal solution Decision variables Objective functions 



n r [rpm] 

P m [kPa] 

s [mm] 

Th [K] 

Tl [K] 

a r H [K] 

A T l [K] 

Nr 

L [mm] 

D c [mm] 

Dr [mm] 

Power [kW] 

VSE [%] 

EA Pi [kPa] 

TOPSIS 

2120 

2550.3 

60.5 

989.6 

298.4 

74.4 

11.8 

339 

70.0 

101.6 

59.5 

6.076 

14.56 

19.690 

LINMAP 

2120 

2550.3 

60.5 

989.6 

298.4 

74.4 

11.8 

339 

70.0 

101.6 

59.5 

6.076 

14.56 

19.690 

Bellman-Zadeh 

2056 

2437.0 

60.5 

989.3 

299.5 

76.4 

12.1 

338 

71.8 

106.1 

58.9 

5.840 

14.51 

18.829 

Experimental 

2500 

4130 

100 

977 

288 

— 

— 

308 

22.6 

22.6 

69.9 

3.96 

12.7 

— 


results [46] 
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Table 3a 

Row data for the sensitivity analysis decision variables to variation of design parameter. 
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Parameter Proposed Variations range Decision variables 






n r [rpm] 

Pm [kPa] 

s [mm] 

Th [K] 

Tl [K] 

A Th [K] 

A T l [I<] 

Nr 

L [mm] 

D c [mm] 

D r [mm] 

A 

1.2 

Lower 

1.2 

2120 

2550.3 

60.5 

989.6 

298.4 

74.4 

11.8 

339 

70.0 

101.6 

59.5 



Upper 

1.5 

2769 

1778.5 

60.7 

981.3 

289.9 

67.2 

6.9 

297 

67.8 

83.5 

56.1 

d(m) 

4 x 10~ 5 

Lower 

3 x 10" 5 

2193 

2485.7 

60.1 

991.2 

296.6 

73.3 

9.1 

306 

72.6 

101.5 

59.9 



Upper 

8 x 10" 5 

1994 

2210.3 

60.4 

984.9 

292.3 

67.5 

7.1 

278 

71.8 

93.6 

59.3 


Note: Sensitivity analysis performed according to TOPSIS method. 


Table 3b 

Row data for the sensitivity analysis of objectives to variation of design parameter. 


Parameters 

Proposed 

value 

Variations range 

Objective functions 

Power [kW] rj S E [%] 

2 A Pi [kPa] 

A 

1.2 

Lower 1.2 

6075.8 

14.56 

19.689 



Upper 1.5 

20897.9 

26.86 

27.053 

d (m) 

4 x 10~ 5 

Lower 3 x 10 -5 

6287.3 

14.39 

19.838 



Upper 8 x 10 -5 

5842.2 

14.74 

17.710 


Note: Sensitivity analysis performed according to TOPSIS method. 


Greek letter 

X ratio of volume during the regenerative processes 

(volumetric ratio) 

er effectiveness of the regenerator 

ri efficiency 

7 specific heat ratio 

H fuzzy membership function 

e emissivity factor 

5 the Stefan’s constant 

v viscosity of the working gas (m 2 s _1 ) 


7. Conclusions 

With regard to previous investigations on Stirling cycle and 
Stirling engine, finite speed thermodynamic has more acceptable 
results in comparison to the others due to considering external ir¬ 
reversibilities. The effects of heat source temperature, temperature 
difference between heat source and working fluid on the heat side, 
engine rotation speed, mean effective pressure, stroke and etc on 
output power, thermal efficiency and system pressure losses are 
scrutinized utilizing NSGAII algorithm. The consequences can be 
implemented for designing and surveying Stirling engine perfor¬ 
mance and robustness. 


Nomenclature 

b coefficient e (0,2) 

C v constant volume specific heat 

Cli decision index in TOPSIS method 

d deviation index 

di + deviation or distance of ith solution from ideal solution 

d z _ deviation or distance of ith solution from non-ideal 

solution 

D diameter (m) 

d wire diameter (m) 

f coefficient related to the friction contribution e (0,1) 

h heat transfer coefficient (W rrT 2 K -1 ) 

m mass of the gas 

N number of gauzes of the matrix, number of regenerators 

per cylinder 
n r rotation speed 

p pressure 

Ap pressure loss 

P output power 

Q. heat transfer rate 

R gas constant 

s stroke 

T temperature 

AT temperature difference 

i ith objective 

j jth solution 

X vector of decision variables, regenerative losses 

coefficient 


Subscripts 

c cylinder, related to the Carnot cycle 

f friction 

H heat source 

L heat sink 

SE Stirling engine 

II related to the second law 

1-4 the processes states 

8 gas 

aver average 

m mean, the system 

app collector aperture 

thrott throttling 
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